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We study the lifetime of the excited state of an atom or molecule near a plane 
surface with a given random surface roughness. In particular, we discuss the impact 
of the scattering of surface modes within the rough surface. Our study is completed 
by considering the lateral correlation length of the decay rate and the variance dis- 
cussing its relation to the Co correlation. 

PACS numbers: 34.35. +a,32.50.+d,68.49.-h,73.20.Mf 



I. INTRODUCTION 

The spontaneous decay rate of an excited atom or molecule is known to depend on its 
environment. This effect is similar to Purcell's effect which has been studied theoretically and 
experimentally in numerous works since the pioneering works by Purcell and Drexhage jll,!^- 
In the very close proximity of a surface the decay rate increases drastically, since the excited 
atom or molecule can couple to non-radiative modes. This can be related to the increase of 
the local density of states (LDOS) near a surface which is due to evanescent modes providing 
more channels into which the excited atom or molecule can decay 

Recently, this effect which allows for controlling the decay rate of atoms and molecules 
has been intensively investigated for random or disordered media. In such materials the 
multiple scattering of electromagnetic modes results in the formation of speckles, i.e., spatial 
fluctuations of the LDOS. Then the spontaneous decay rate of atoms or molecules close or 
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within such systems becomes a statistical quantity, which depends on the one hand on 
the local near-field environment of the source and on the other hand on the mesoscopic 
fluctuations of the random material itself. In particular, the fluorescence rate statistics 

nn 

or fluctuations of the LDOS in such media has been considered theoretically |5H2[ and 



experimentally 
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A random rough surface is similar to a bulk disordered medium in the sense that above 
such a surface the LDOS shows a spatial spech.e pattern Q. The Metitue of au atont or 
molecule becomes a random quantity which depends on the local environment of the particle 
and the statistical properties of the surface. Recently, the speckle pattern above random 



media 



15Ml7j| has been studied. The goal of this work is to reconsider the impact of surface 



roughness on the spontaneous decay rate. Previous studies have considered the impact of the 
surface roughness on the average decay rate for atoms or molecules near metal surfaces [lsL 



2l| . Here, we will focus on the lateral correlation of the decay rates, i.e., the correlation 



between the decay rates of an atom or molecule placed at different positions above the rough 
surface by keeping the distance to the mean surface constant, and its variance. In addition, 
we will specifically consider the influence of surface modes for which the enhancement of the 
decay rate is very large. 

For pedagogical reasons we consider a semi-infinite SiC material with a rough surface. 
Firstly, we do not have to consider nonlocal effects for the considered atom-surface distances 
which can be quite important for metal surfaces 



22|-|25|. Secondly, SiC has only one well 

for its permittiv- 



established surface resonance and is well described by a simple model 
ity Nonetheless, the results we obtain are applicable to arbitrary local, homogeneous and 
isotropic materials. For the description of the surface roughness we use the perturbation 



theory introduced in Ref. 



27| up to second order within the surface profile function. 



The paper is organized as follows: In Sec. 2 we give a short introduction to the calculation 
of decay rates close to a plane surface in the weak coupling limit. The effect of the surface 
roughness on the mean decay rate is discussed in Sec. 3 where we also give some interpretation 
of the observed roughness correction. Finally in Sec. 4 we investigate the lateral correlation 
and the fluctuations of the decay rate. We finish with the conclusion in Sec. 5. 
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II. SPONTANEOUS DECAY RATE 
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For an electric-dipole transition in the weak-coupling regime, the normalized spontaneous 
decay rate of an atom or molecule placed at r A can be expressed as [28] 

li 

where ej is the unit vector in the direction of the dipole transition, t symbolizes the trans- 
posed vector, fc = u /c with u the frequency of the dipole transition. is the decay rate 
in free space (see for example in Ref. 28|). G is the classical electric Green's dyadic for the 
geometry considered. In our case, this geometry consists of a half-space of a given material 
characterized by its permittivity e(u) with a rough surface as depicted in Fig. (TT]). 




Figure 1: Sketch of a dipole at r A in a distance z A above the mean of a rough surface described 
by a profile function S. 

Before considering the role of surface roughness, we summarize the known results for 
a flat surface. We can derive the decay rate from Eq. ([1]) for a dipole moment parallel 
Tn and perpendicular T± to the surface by inserting the Green's dyadic from Eq. flAll) in 
appendix [A] We find the well-known relations 

rf 3 r k °d> 



r 4 



Here, we have introduced the lateral wave vector n = (k x ,k y y, the perpendicular wave 
vector 7 r = i/fc^ — k 2 and 7 = ^ k 2 — k^; r s and r p are Fresnel's coefficient for s- and p- 
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polarized light. Additionally, we have split the decay rate into its radiative {k < ko) and 
non-radiative (« > k ) part. 
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Figure 2: (Color online) Normalized decay rate of an atom in a distance za above a flat surface 
with its dipole moment oriented (left) parallel and (right) perpendicular to the surface. We have 
chosen the material properties of SiC at the transition frequencies ujq = 1.4- 10 14 s -1 , 1.787-10 s -1 
and 3- 10 14 s" 1 . 



In Fig. El we show a plot of the distance dependence of the decay rates F\\ and T± for 
an atom near a flat SiC interface. We consider a transition frequency coinciding with the 
SiC surface phonon resonance at Waj SPhP = 1-787- 10 14 s _1 and a transition frequency slightly 
smaller and slightly larger than wsphP- In all three cases one finds the known characteristics 
of the decay rate near a flat surface [4|: i) For relative large distances za > \ = 2txc/ujq the 
decay rate oscillates due to the phase change of the reflected field, ii) In the near-field regime 
with za < Ao = 2ttc/uo the decay rate is highly increased due to the decay into nonradiative 
or evanescent channels. For u — 1.4 • 10 14 s -1 and u = 3 • 10 14 s _1 the atom or molecule can 
decay into total internal reflection modes, whereas for u = 1.787 • 10 14 s _1 it can decay into 
surface phonon polaritons. It can be expected that within this near-field regime, the decay 
rate will be very sensitive to the multiple scattering of surface waves within a rough surface 
for oj = 1.787 • 10 14 s" 1 . iii) Finally, for a distance za smaller than about 100 nm the decay 
rate diverg . This so-called quenching effect emerges from the 1/z 3 electrostatic 

interaction of the atom's dipole field with the surface. Therefore, this effect is extremely 
localized so that in this extreme near-field regime the decay rate will only be sensitive to the 
change of the local environment of the atom or molecule as for example to the local change 
of the surface geometry due to roughness. 
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III. ROUGHNESS CORRECTION TO THE DECAY RATE 

Now, we turn to the effect of surface roughness on the decay rate. To this end, we 
consider a stochastic surface profile function 5" describing the deviation of the rough surface 
from flatness (see Fig. [1]). The function S is modeled as a stochastic Gaussian process with 
mean value and correlation function given by 

<S(x)> = 0, (4) 

(S(x)S(x!)) = 5 2 e~^^ = 5 2 W(\x - x'|), (5) 

x = (x,yy. The brackets () stand for the average over an ensemble of realizations of the 
surface profile S(x); 5 is the rms height and a the correlation length of the surface profile. It 
follows that the Fourier components S(k) of the surface profile function fulfill the relations 

(§(*)) = 0, (6) 
(S(k)S(k')) = (27r) 2 5 2 5(n + ( 7 ) 
where we have introduced the surface roughness power spectrum 

g ( K ) = Jd 2 xW(\x\)e~ iK - x . (8) 
In the following calculations we will assume a gaussian correlation function W(|x — x'|) = 

K?a 2 

exp(|x — x'| 2 /a 2 ) so that in this case g(n) = 7ra 2 e ~ . By introducing a stochastic surface 
profile, the fields are scattered by that surface and hence the decay rate becomes a stochastic 
processes. The reflected fields can be described by a stochastic reflection coefficient which 
determines the decay rate in Eqs. (J2J) and ([3]). The statistics of the decay rate is itself 
determined by its mean value and higher moments. Here, we will concentrate on the mean 
decay rate and in the next section we will turn to the correlation function and the variance. 
By virtue of Eqs. (J2J) and ([3]) the mean decay rate depends on the mean reflection coefficients 
of the surface. The average restores the translational invariance so that the decay rate 
depends only on the distance to the surface za- The mean reflection coefficients can be 
determined perturbatively if the surface roughness is much smaller than the wavelength 



Ao- It has been shown in Ref. 29j that by using the perturbation theory of Ref. 27|, the 
correction to the Fresnel reflection coefficient Ar s / P = (r s / p ) — r s / p due to roughness is up to 
second-order in the surface profile given as 

Ar s/P = -2i 7r (A/ P ) 2 M s/p (9) 
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where 



7r + 7t 



and D r 



lr 



1^ + 7 t 



(10) 



using 7t = \jk\e — k 2 . The expressions for M s and M p can be found in Ref. 29]. Therefore, 
one can easily get the second-order correction to the decay rates 
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by replacing the reflection coefficient in Eqs. ([2]) and fl3]) by Ar s and Arv Now, we use the 



approximation for the correction to the reflection coefficient from Ref. 

Ar s/p « r s/p 2K 2 5 2 



29] 



(12) 



which holds in the quasi-static regime for k ^> ko and /ta > 1. Since this is in the quasi-static 
regime equivalent to distances za <C a, we can conclude that for za <C a we have 
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(13) 
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Figure 3: (Color online) Correction to the decay rate Ar for an atom with its dipole moment 
oriented (left) parallel and (right) perpendicular to the mean surface. The parameters are the 
same as in Fig. [2j For the rough surface we choose 5 = 5nm and a = 200 nm. 

In order to illustrate the effect of roughness we plot in Fig. [3] the roughness correction 
ATy/i. for the same frequencies as in Fig. [2] considering a rough surface with an rms 5 = 5 nm 
and a correlation length a = 200 nm. It can be seen that the roughness correction is 
very small in the large distance regime for za > Ao, but can be relatively large for small 
distances, i.e., for za < A where the decay rate is very large due to the decay into non- 
radiative channels. As will be discussed in more detail in the following, electrostatic effects 



7 



in the extreme near-field for za < 100 nm and surface phonon polaritons in the intermediate 
distance regime are responsible for this relatively large correction. The first is a local effect, 
whereas the latter is a multiple scattering effect. 

In the intermediate distance regime 100 — lOOOnm it can be seen on Fig. [3] that the 
roughness correction is slightly positive at a distance of ~ 1000 nm when total internal 
reflection modes are excited (i.e., for lo$ = 1.4- 10 14 s _1 and Uq = 3- 10 14 s -1 ). For a frequency 
ujq = wsPhP surface phonon polaritons are excited in this distance regime. Surprisingly, 
the presence of roughness leads to a large negative correction indicating that the lifetime is 



increased. This effect has been studied in Ref. 



29| in terms of the LDOS (see Fig. 14 Ref. [29|) 



and is due to the roughness induced multiple scattering of surface modes. The scattering 
causes a broadening of the dispersion relation 311 ] . Due to this broadening as illustrated in 
Fig. IH the LDOS (which is proportional to Im(r p )) becomes smaller for frequencies close to 
c^sphP for intermediate k or distances z&, respectively, explaining the observed decrease of 
the decay rate. 
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Figure 4: (Color online) Plot of Im(r p ) for the flat and Im((r p )) for the rough surface for k 
3.3 • 10 6 m _1 which corresponds approximately to a distance of za ~ K ~ l = 3 • 10~ 7 m . 



We now consider the very small distances za a regime. The roughness correction is in 
this case due to the local electrostatic interaction (quenching) of the atom dipole moment 
with the rough surface resulting in a positive and large roughness correction to the decay 
rate. This quenching effect is similar for all transition frequencies as can be seen in Fig. |3J 
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This correction can be described by the quasi-static expression in Eq. (lT3"j) . We will now see 
that we can retrieve this expression in Eq. f[T3"j) using a simple physical argument. If <C a, 
curvature effects are negligible so that the atom feels only the local deviation of the surface 
from flatness. This effect can be described by the ansatz: 



(rn /x )«(r(d + 5(x))>. 



(14) 



This means that one replaces locally the surface profile by a shifted flat surface. Employing 
this approximation in Eqs. (T5]) and ([3]) we see that it is equivalent to the replacement of the 
mean reflection coefficient by 

(r s/p )=r s/p (e 2i > 5 W). (15) 

Note, that this is the expression for the propagating (n < ko) and the evanescent (k > ko) 
part. Now, this expression can be easily evaluated, since we have assumed that S(x) is 
Gaussian distributed. We find 

(r s/p )=r s/p e- 2 ^ 52 . (16) 
For propagating waves such that k < ko, this is the well-known result of the Kirchhoff 



approximation 



30| which holds if A <C a. But for k ^> ko, this approximation produces the 



result in Eq. (112]) when expanding the exponential up to second order in the rms S. Hence, 
we have retrieved Eq. (ITS]) . 



IV. STATISTICAL PROPERTIES OF SPONTANEOUS EMISSION 



We now examine the fluctuations and the spatial correlations of the decay rate above 
a random rough surface. The statistical properties of fields in random media has received 
a lot of attention in the past thirty years. Here, we are interested in some recent results 
relevant for our system. It has been shown recently that the intensity correlations above 
random media or materials with rough surfaces become non universal in the near field, i.e., 
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-|34|. 



they highly depend on the properties of the random media or rough surfaces 
A remarkable connection has been established between the LDOS fluctuations and the Co 
correlation js, 7| for multiple scattering media. The Cq correlation is defined as the infinite 



range contribution to the correlation of the intensity in multiple scattering media [6j. A 
simple explanation has been reported recently Finally, the multiple scattering of surface 



modes in a random media or rough surface can lead to localized surface modes 



35 



which 



9 



show a characteristic long tail distribution of the intensity enhancement of the fields close 



to the surface 
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This effect can be neglected for surfaces with small roughnesses as 



shown for fractal surfaces [39|, (^(J. Here, we focus on the lateral correlation and the variance 
of the decay rate above a rough surface and we will discuss the relation of the variance to 
the LDOS fluctuations and Co correlation. Since the perturbative approach is restricted to 
small surface roughnesses, we will leave the problem of localization and its relation to the 
distribution of decay rates for future studies. 



A. Variance and correlation function 



Before evaluating the correlation function of the decay rate, we first determine the vari- 
ance which is up to second order in the surface profile given by 

= (it) - <r 4 > 2 = <rf rf >. (17) 

Obviously the variance is a special case of the more general correlation function 

(r i (r)r J -(r')) = (r.r^pee + (r.r;.)^, (is) 

which can be divided into a specular (depending only on the mean field) and a diffuse con- 
tribution (due to the fluctuating part of the field). By inserting the perturbation expansion 
Tj « 1^ + + we find up to second order for both of these contributions 

(r,r;.) spec = rf^rfV) + rS 0) (*)<rf V)) 

(19) 

+ <rf(z)>rf (zf) 

and 

(r i rj>d lff = <r, (1) (r)rW(r')>. (20) 

Obviously, the specular part only depends on z and z' . This is due to the fact that for 
the mean field the translational symmetry with respect to the x-y plane is restored after 
averaging, whereas for a fixed z and z' the diffuse part contains lateral correlations with 
respect to |x — x'|. Furthermore, we note that the variance depends on the diffuse part of 
the correlation function, only. 
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B. Lateral Correlation 

Now, we focus on the lateral correlation only. Therefore, we assume z = z' = za- The 
correlation function is then given by 

(r i (x)r i (xO) = (r,r,) spec + (r,(x)r i (x')) dJff . (21) 

It is clear from Eq. (fl9|) that the specular contribution is just a constant term giving an 
infinite range correlation depending on the distance z\ only. On the other hand, the diffuse 
part depends on |x — x'| and therefore contains the lateral short range correlations. We 
derive the explicit expression for the diffuse part of the correlation function in appendix O 
The result can be stated as [see Eq. (jC6D ] 

(rS^gV)) = (|f 2Re J^ gmm . Za)gm . 2A)ei , (x -xo (22) 



where the functions Gi and Fj are defined in Eq. flC7j) and (1C8|) . In the following we will 
discuss this expression in more detail. 

Let us focus on the evanescent regime, i.e., z^ <C A. Then the exponential function 
exp(i7^£A.) ~ exp(— \k ± £\z&) in the integrand <Zj of Fj and Gi [see Eq. (102[) ] acts as a low 
pass filter and restricts the contributing £ to £ < 1/za- On the other hand, for a Gaussian 
roughness correlation the roughness power spectrum g(|£|) also acts as a low pass filter 
restricting the £ to £ < 1/a. Therefore we can make simple approximations for Eq. (1221) in 
the two limits a z^ and a ^ z\. 

In the case a ^> the functions Fj and Gi can be approximated by Fj(0;zA) and 
Gj(0; Za)- It follows immediately from Eq. (jC6ft 



/ r (l) r (l)'\ r 2 

V V ~ (37r) 2 4Im[G i (0; ^a)]^^^; ^a)]t2^(|x - x'|). (23) 

1 oo ft 

We can conclude from this expression that the lateral correlation function reproduces the 
correlation function of the surface roughness for distances such that z& a, which also 
holds for za > A. In particular the correlation length of the lifetime correlations coincides 
with a. This means that in this distance regime one can directly measure the correlation 
of surface roughness by measuring the correlations of lifetimes above the surface. In the 
quasi-static limit the correlation function can be further simplified (see appendix [D] for a 
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detailed calculation). We find 



(r 



ii/jJii/jjdiff 



r (0) (z*)r {0) 



5 2 

9° W(\x-x! 



(24) 



Now, in the opposite limit, where we have a < za the roughness power spectrum g(|£|) 
in Eq. (jHJ) can be approximated by g(0) = na 2 and can be taken out of the integral. The 
remaining expression can be further simplified in the evanescent regime assuming that the 
most important contributions stem from £ 3> ko and £ ^> ko\e\. The resulting expression for 
the r || can be written as (see Eq. (1D17[) appendix [Dj) 



3 V V Z A + (I X - X 'I 

XV(r)ry(r))oc5 a + 



/I212 ' 



(25) 



where P3 is the Legendre polynomial of third power. Hence, for distances such that a <C 
*C A the lateral correlation goes rapidly to zero for |x — x'| ^> za- Surprisingly, the 
lateral correlation length only depends on z^ and does neither depend on the correlation 
length of the surface roughness nor on the properties of the material. For ^r^ 1 '* (x)r^ 1 ' 1 (x')) 
we find a somewhat more complicated but similar expression in Eq. flD18j) which leads to 
the same conclusions. We note that a similar result was found for the intensity correlation 



in the near-field of a random medium 
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Figure 5: (Color online) Plot of the correlation functions (rj^(r)rY^(r')^ and <^rj'| ±y ' (r )rj~| ±yl (r') ) in 
Eq. (|25p and (|D18j) over the lateral distance (x — x')/za normalized to their value at x = x'. 
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In Fig. O we plot the quasistatic results of ^T^ L ' ) (x)r^(x / )) and (T^P (x)rf| (x')) from 
Eq. (ID17P and (1D18|) for a fixed distance za which is assumed to be so small that the 
conditions for the quasistatic approximations are met, but za S> cl- Note that the regime 
where these approximations are valid might be hard to achieve in practice, since the non- 
retarded regime starts for SiC for example for za < 200 nm, i.e., for distances which are 
not much larger than typical surface roughness correlation lengths. Nonetheless, one can 
expect that, for intermediate region, the correlation length is the atom-surface distance z\. 
To illustrate this fact, we plot in Fig. [H] numerical results for the correlation function in 
Eq. (|2~3j) . It can be seen that for intermediate distances a < za < A the correlation length is 
about za- 
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Figure 6: (Color online) Plot of the correlation functions (left) (r^(r)r^(r')) and (right) 
(r|, 1 ^(r)r|| 1 ^(r / )) in Eq. (j23|) for different distances z at ujq = 1.787 • 10 14 s~ 1 choosing the surface 
roughness parameters 5 = 5 nm and a = 50 nm (|D18j) . The correlation functions are normalized to 
their value at x = x'. Furthermore we plot again the approximations shown in Fig. [5j 



Variance and standard deviation 



Let us now return to the variance given by Eq. 017p . In the quasistatic limit using Eq. 
we find for z <C a the simple expression 



IIA 



r (0) (z) 
1 haw 



-i 2 



(26) 



In Fig. [7] we show some plots of the standard deviation cr/T . It can be seen that the 
standard deviation approaches the approximate result 35 /za for za <C a so that for very 
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small distances as z = 10 nm the standard deviation or variance is on the order of 1^°) 
or resp., indicating that fluctuations of the decay rate are large in the quasi-static 

regime. In all shown cases a falls off rapidly with the surface atom distance, which is due 
to the small roughness considered. It could be expected from Eq. (fSB]) that in the distance 
regime a C 2a the standard deviation varies like a oc 5a/ ' z\. Indeed, in Fig. [7] we find this 
power law for a± but not for au. Note, that for the distance region around 200 nm for which 
we have found a relatively large roughness correction of 15% to the mean decay rate cr±/\\ is 
smaller than 10%. 
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Figure 7: (Color online) Plot of the standard deviation a±_ (left) and cry (right) for different distances 
of an atom placed at the distance z\ above a rough SiC surface with 5 = 5nm and a = 200 nm. 
The transition frequencies are chosen to be ujq = 1 • 10 14 s , 1.787 • 10 14 s _1 and 3 • 10 14 s . 



D. LDOS fluctuations and Co correlation 

Finally we want to explore the relation between the LDOS and the infinite range intensity 



correlation Co as studied for multiple scattering media It was shown that this infinite 
range correlation equals the LDOS fluctuations jjj]. Since the decay rate is proportional to 
the LDOS, one could expect by analogy that above a rough surface the C correlation equals 
the decay rate fluctuations or more precisely Cq = <j 2 /(T) 2 . To prove this, we follow the 



reasoning of Ref. js] adapted to our problem. First we assume for simplicity that we have 
a non-absorbing half-space with a rough surface. Then the radiated power of a dipole in 
vacuum Pq and of a dipole above the rough surface P can be related to the decay rate of an 
atom in vacuum T and the decay rate of an atom above a rough surface T by the simple 
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relation 28 1 

p r 

f = jr ■ (27) 
Since the decay rate is proportional to the LDOS p at the position of the atom we also have 
P/P = p/p . 

Now, we define the speckle correlation function C(x, x') for the lateral correlations with 
respect to x and x' above the surface as 

_ (/(x)J(xQ) _ 
C(X ' X) - (J(x))(/(xO) ( } 

where /(x) is the radiated power of the dipole at the position x such that the integral over 
a plane parallel to the mean surface at z = (but for z > z& and z > \ so that evanescent 
waves are not included) 

lim \ I d 2 x/(x) = P-R (29) 

r ->oo A J A 

gives the total radiated power into the halfspace for z > za- Here, we have introduced the 
circular area A = irr" 2 . The factor R takes into account that only a part of the total power 
P is radiated into the halfspace for z > z\, whereas another part P ■ T is radiated into the 
halfspace z < such that the total power is the sum of both contributions, i.e., T + R = 1. 
With the relation P/Pq = p/po and the definition of /(x) and C(x, x') we have 

y Q / ro^oo A 2 R? Pq 



r 2 [d 2 x [ dV (/(X)J(X')) 
J A J A 



(30) 



'■()- 

Since after averaging we retrieve the translational invariance parallel to the plane with z = 
and isotropy, we have (/(x)) = P ■ R = P - Rp/p . Inserting this relation into (15U|) it follows 
that 

(p 2 ) = (p) 2 lim -L f d 2 x / dV[l + C(|x-x'|)] 

ro— >oo /{ J A J A 

\ 2tx f°° 1 

= (p) 2 1+ lim — / drrC(r) 

ro->oo A J 

In fact, the integral over C(r) gives nonzero contributions only for C(r) = const = Co, 
when assuming that for intensity correlations the relation C(0) > C(r) for r > is valid. 
Therefore the integral in Eq. (T5TT) reduces to Co, i.e., the constant component of C(r) which 
is the searched for infinite range Co correlation. Hence, we find 

r _ (p 2 ) - (p) 2 _ (r 2 ) - (r) 2 
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which proves our statement that the Cq correlation equals the decay rate fluctuations or 

c = <x 2 /(r) 2 . 

Finally, using the perturbation expansion for the decay rate, we get for the Co correlations 
above a rough surface 



a 



2 



(rWrW) (rWrW) 1 (rWrW) 



Co <r> 2 ^ (r(°) + (r(2)» 2 r(°) 2 / r(°) 2 " (33) 



1 + r (o) 

Hence, the square of the normalized standard deviation a/T^ gives the infinite range Co 
correlation showing its sensitivity to the local environment which enters through its depen- 
dence on Za- In the quasi-static regime, we can now use the above derived result in Eq. 
so that for C o we find Co ~ 96 2 / z\. We note that this result is very similar to the result 



found by Shapiro js] for random media, where Co = ir/(kl) with the wave number k and the 
mean free path of radiation inside the random medium /. In our case, S 2 /z\ corresponds to 
the scattering strength l/(kl). 



V. CONCLUSION 



We have studied the impact of surface roughness on the decay rate or inverse lifetime of 
a molecule or atom above a rough surface. For pedagogical reasons we have only considered 
SiC as the bulk material, but the conclusions can be easily transfered to other dielectric 
materials supporting surface modes as for example silica. Our results show that the decay 
rate might be reduced by 15% due to the surface roughness for very shallow roughnesses 
with an rms of 5 = 5 nm and a correlation length of a = 200 nm. This reduction is due to 
the surface induced scattering of surface modes which prevails for intermediate distances. 
On the other hand, for very small distances za <C a the rouhgness correction to the decay 
rate r'°) is due to the local electrostatic interaction of the atom or molecule with the surface 
and is given by the simple expression 65 2 / z A T^°\ 

In addition, we have studied the variance and the lateral correlation of decay rates above 
the rough surface. We find that the lateral correlation length is approximately given by the 
distance za itself in the nonretarded regime for distances larger than a. For distances smaller 
than a the lateral correlation function resembles the surface roughness correlation function 
allowing for a direct measurement of the surface roughness properties by measuring decay 
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rate or lifetime correlations. The variance itself is a special case of the lateral correlation 
function and we have pointed out that it equals the Co correlation as for random media. 
We have shown that it can also be approximated by a simple result o 2 = 95 2 /z A in the 
quasistatic regime for za <C a emphasizing that the infinite range Co correlation highly 
depends on the local environment of the atom, i.e., on the distance za- 
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Appendix A: Green's function for a fiat surface 



The Green's function with observation point and source point above the flat surface, i.e., 



for < z < z' can be stated as 
G (0) (r,r» = 



11 



43] 

d 2 K ie iK -( x " x ') 
1 

3k 



t__ e ^{z'-z) + to e M*'+*) 



(Al) 



^S(z - z')S(x - x')e 2 g> e 2 



'0 



where e z is the unit vector in z-direction and ® symbolizes the dyadic product. The tensors 
1 and R are defined as 



1. 



K K) &j [K) 



j={s,p} 
j={s,p} 



(A2) 
(A3) 



where 



a s ( K ) = a+O) = -(-k v , k x , 0)\ 

Hi 



a p [K 



2\t 



(A4) 
(A5) 
(A6) 



are the polarization vectors for s- and p-polarization. Note that these vectors are always or- 
thogonal, but only normalized for propagating modes with k < k . The reflection coefficients 
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r s and r p are the usual Fresnel coefficients 

r s = and r p = . (A7) 

7r + 7t 7r£ + 7t 

Appendix B: First-order Green's function 

The correction to the Green's function we find from first-order perturbation theory is [2 

G«(r,r» = - f^L fc °' (e -V §W( K > - K)e^'^X + _(K, *') (Bl) 

J {2iry J [2iry 4-f T % 

with 

S^(k' -k) = Jd 2 xe- i{K '~ K) *S(x). (B2) 

and 

X+_(#e,K') = k i( K> > ® »i ( B3 ) 

*»i={s,p} 

The elements of the tensor are given as 

X ss = i^/s ■ k' (B4) 

X sp = -t s t p -^-e 2 ■ k X k' (B5) 
V efc o 

^ P s = -t p t' s -^—e z ■ kx k' (B6) 
V efc o 

^pp = +tptJ,-T2(K«'e - 7t7t* ■ k ') ( B7 ) 
e/c 

where k— k/k and t s , i p are the usual amplitude transmission coefficients 

t, = ^- and t, = (B8) 
7r + 7t 7 r e + 7t 

Appendix C: Correlation function 

By inserting the Green's function from Eq. (IBlj) into Eq. (JT]) we find for the first-order 
correction to the decay rate the expression 

(CI) 



Too ^0 

where 



^o(e Jlgi^+^^.x^-eJ. (C2) 
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With this definition at hand it is an easy task to check that the correlation function is 
(rW(r)rf (r')> ( 37 r)2 f d 2 „ f dV 



r 2 P 

1 oo ft 



using the relations 



K + = K " + (k-k'), (C4) 
«- = «"-(«-«')■ ( C5 ) 



By introducing the new variable ^ = k - k' we can write the correlation function as 



r 2 P 

1 oo ft 



2Re j^&gmmrtGfa*)^-^ (C6) 



with 



and 



d 2 /t 

= I j2^ a i( K , K - & z ) ( C7 ) 



F^zO = (C8) 



Appendix D: Approximations for quasi-static limit 

In the quasi-static limit (« ^> ko) the reflection coefficients can be approximated by 

r p Rj -— and r s « ( D1 ) 

By inserting these relations into Eqs. (IC7I) we find the quasi-static approximations for the 
decay rates 

rf|0) 3 m 



^ 16 (k z) 3 \e + l 
r ( , 0) 3 1 . /e-1 



Too 8(M 3 Ve+1 



Im . (D3) 



In particular, = 2T^\ 
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i) distance regime z <C a 



Now, we want to find similar simple approximate expressions for the correlation function 
in Eq. (1C6|) . To this end we consider first /t>(, which is fulfilled for z <C a. For such wave 
vectors we can approximate Eq. flC7j) by 



dK feg(e-l) 2i7iZ 
2tt 872 



t 2 - 2 
} 

6/Cn 



t 2 _ % 



[K 2 (e+l)-kle] 
l)-k 2 e]. 



2tx 4 7r 2 

Using the quasi-static approximation for the transmission coefficients 



to ~ 1 and t T 



e+ 1 



allows for further simplification. We find 

G,|(fcz) w 
G x 



3! 1 1 e - 1 
'4^(2^) 4 fcfe + 1 
3! 1 1 e - 1 



2tt (2z) 4 /c 2 e + 1 
Inserting these approximations into Eq. ( 1C6|) finally yields 

<r[ | 1) (r)r[ l 1 V')> si 5 2 r / 



r 2 



256 k%z^z' A 



Im 



Ir^wr^r')) 81 5 2 



r 2 



64 fc^V 4 



Im 



e - 1 



e- 1 
e+1 



^(|x-x'|), 
^(|x-x'|). 



As can be expected from G± = 2G» we find 



<r«(r)r«(r')> 



(rj 1) (r)rj 1) (r')) 



(D4) 
(D5) 

(D6) 

(D7) 
(D8) 



(D9) 
(D10) 



(Dll) 



ii) distance regime z a 



In this limit, we consider the case k ^> £ yielding 

d 2 K 



G11 



(27T) 2 



(D12) 

(D13) 
(D14) 
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Together with the quasi-static approximation, i.e, £ ^> k and £ ^> k e we get 

G± » e" ?2 £ 



^::V .M E/l(z) e^ 



2tt /co4i7 r 



e + 1 



(D15) 



and 



1 f dK A; 2 (e - 1) 



Gii ^e-^cos 2 (^)- / — k 



272,- 4i 7r - L *jj e *gv£J (D16) 

where we have introduced £ x = £ cos and £ y — £ sin(6 l ). Finally, when plugging these results 
into Eq. we find [42] 



(r^(x)rVV )) 
r 2 , 

(rj^rfiV)} 5 v 



'1212' 



x — X 



2-^1 ( 2, 2 > 3; 



+ (x - x') 2 
Ix-x'l 2 



(D17) 



20 |x - x 

T — 



'12 



2-M ( 3, 2' ^' 



x — X 



'|2 



(D18) 



where P3 is the Legendre polynomial of third power and 2 F x is the hypergeometric function. 
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